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M/G/l  with  Exceptional  Service  and  Arrival  RaU 


Martin  Krakowski 


Abstract 


-^The  model  M/G/l  is  modified  by  to  providing  the  pioneer  customer,  i.e.  the  customer  who 
terminates  an  idle  period  (and  initiates  a  busy  period)  with  exceptional  service;  and  (Q)  by  allowing  an 
exceptional  arrival  rate  during  the  idle  period.  The  regimen  is  FCFS  and  the  server  idles  only  when 

customers  are  absent.  Ml*  derived omni-equations  for  the  delay,  for  the  backlog,  and  for  the  queue  size 

f  f r,  ■  ■  ’ 

as  found  by  real  or  virtual  arrivals.  relate; these  processes  to  the  regular  M/G/l  as  convolutions  of 
the  delay  in  M/G/l  with  modifying  variables  in  the  model  treated.  The  queue  size  is  derived  from  the 
delay  by  applying  the  Poisson  operator. 

This  paper  extends  the  results  of  P.D  Welch  (1964)  who  allowed  exceptional  service  but  not 
exceptional  arrival  rate,  and  who  did  not  discuss  composition  relations.  His  solutions  used  Laplace 
transforms  and  generating  functions.  Wo  -believe  that  the  current  treatment  is  simpler  and  more 


suitable  for  further  generalizations.  Introducing  the  exceptional  arrival  rate  forces  us  to  distinguish 


between  the  perception  of  the  observer  and  of  the  customer  source. 
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Notation 


<r=  arrival  rate  when  server  idles 
A  =  arrival  rate  when  server  serves 

y  =  service  time  for  a  “pioneer"  (customer  who  initiates  a  busy  period ) 
x=  service  time  for  non-pioneers,  or  “followers" 
z=residual  service  time  of  z  or,  briefly,  the  residue  of  z 
p  —  AEx  =  Ap  Po~ff  y 

B=backlog  or  unfinished  work 

w  =  overall  delay  as  perceived  by  the  customer  source 

w„=  positive  delay,  i  e.  delay  of  a  follower;  a  pioneer’s  delay  =  0 

u=  clearance  time  of  the  system  (=  virtual  waiting  time  in  our  model) 

u,  =  clearance  time  of  the  system,  provided  the  server  works 

n  =  queue  size  at  a  random  instant  (or  viewed  continuously) 

n*=queue  size  at  a  random  instant  provided  the  server  works 

qd=queue  size  left  behind  by  a  departing  customer  who  enters  the  service  station 

qa=queue  size  found  by  an  arriving  customer 

q<(j  =  queue  size  left  behind  by  a  departing  follower  who  enters  the  queue 
q*a=queue  size  found  by  a  follower  (i.e.found  by  a  customer  while  server  works 
N=system  size  as  seen  by  a  continuous  (or  Poissonian)  observer 
Q  =  system  size  found  by  a  newcomer  (or  left  behind  by  a  served  customer) 

Q0=fraction  of  newcomers  who  find  the  server  idle 

Q*=fraction  of  newcomers  who  find  the  server  busy 

Po  =  Pr(N  =  0)  P„=  Pr(N>0)  =  Pr(server  serves)  =  1—  P0 

A  free  copy  of  a  random  variable  z  is  a  random  variable  having  the  same  distribution  as  z  and 
independent  of  any  other  variables  within  the  same  argument.  Different  free  copies  of  the  same 
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variable  are  usually  designated  by  the  same  symbol  if  they  occur  in  different  arguments:  this  causes  no 
confusion. 

Omni-  Transform 

Definition  1  The  omni-transform  of  the  random  variable  A  is  the  expectation  of  an  arbitrary 
function  of  A:  E0( A).  The  arbitrariness  of  ii>  is  limited  only  by  the  requirement  that  Ev>(  A)  and  v;( 

)  exist.  When  using  this  definition  we  usually  apply  the  so-called  omni-convention  (see  below). 

Definition  2  The  omni- transform  of  the  random  variable  A  is  an  arbitrary  functional  of  A. 

It  can  be  shown  that  the  two  definitions  are  essentially  equivalent  since  under  wide  enough  conditions  a 
functional  of  A  can  be  represented  in  the  form  Eip(A.)  for  some  xl>.  Among  the  advantages  of  Definition 
2  is  that  we  need  not  know  the  function  whose  expectation  is  the  functional  of  interest  and  we  can  thus 
easier  specialize  ii>( A).  Thus,  if  we  need  Pr(A<t)  we  need  not  write  v(A)=EH(t  — A)  where 
H(t  — A)=l  if  t<  A  and  =0  otherwise;  we  simply  write  ij>(A)  =  Pr(  A<t).  Another  advantage  of 
Definition  2  is  that  no  need  arises  for  the  omni-convention  (see  below).  But  perhaps  Definition  1  is 
somewhat  simpler  when  balancing  ii>( A).  Since  the  omni-equations  are  typographically  identical  for 
both  definitions  (when  the  omni-convention  is  used  with  Definition  1)  the  reader  may  choose  either. 

Omni- Convention  Omni-equations  are  easily  told  by  sight.  Hence  Et^A)  can  be  replaced  in  print  by 
i/.(A)  and  the  operator  E  is  retained  mentally;  this  is  the  omni-convention.  Definition  2  has  no  need  for 
the  omni-convention.  Omni-equations  look  the  same  under  both  definitions  if  the  omni-convention  goes 
with  Definition  1.  One  can  even  jockey  between  the  two  definitions  for  didactic  or  esthetic  reasons. 
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Section  i.  The  Balance  of  t/>(B)  in  G/G/l  with  Exceptional  Service 

It  is  of  interest  to  see  how  far  we  can  get  within  a  modified  G/G/l,  in  which  the  pioneer  customer 
receives  exceptional  service,  when  we  apply  the  method  of  omni-equations.  (The  generality  of  the 
arrival  stream  includes  the  possibility  of  exceptional  rate  for  the  pioneers.  )  Once  we  reach  a  barrier  to 
further  progress  we  will  see  clearer  why  that  barrier  may  be  lifted  by  a  fundamental  property  of  a 
steady  Poissonian  source,  namely  that  true  arrivals  see  the  same  picture,  stochastically,  as  does  a 
continuous  observer  or  a  Poissonian  observer. 

In  our  model  the  backlog  B  manifests  itself  to  new  arrivals  as  their  delay  w;  we  can  say  that 
“w  =  B  conditioned  upon  a  customer  just  to  arrive”  and  it  manifests  itself  to  a  continuous  observer  as 
the  virtual  delay  u  (or  clearance  time  in  our  model).  An  observer  who  samples  the  system  at  instants 
generated  by  an  independent  Poisson  source  of  steady  intensity  is  equivalent  to  a  continuous  observer. 
Clear  intuitively  that  Poissonian  sampling  should  be  unbiased,  this  was  also  treated  analytically  (  Wolff 
1982).  Such  an  observer  will  be  tersely  referred  to  as  a  random  observer  and  his  observations  as 
random  observations.  We  can  say  that  “u  =  B  conditioned  upon  the  observation  instant  being 
random.”  The  process  B  is  thus  a  superposition  of  the  point  process  w  and  the  continuous  (almost 
everywhere)  process  u . 

Under  steady-state  conditions  B  is  a  balanced  process,  which  means  that  E(dB)=0  during  a 
random  dt.  But  ti’(B),  an  arbitrary  function  of  B,  is  also  balanced  since  we  also  have  dtl>(B)=0 
over  any  random  time  interrval  dt. 

We  shall  now  derive  the  balance  equation  for  cKB).  Balancing  an  arbitrary  function  of  a 

process,  rather  than  the  process  itself,  is  the  essence  of  the  omnt-  method.  The  expected  changes  in 
u>(B)  derive  from  two  causes: 

(a)  While  the  server  works  the  backlog  B  is  being  worked  off  continuously  at  the  rate 


du=du»=— dt;  and  du  =  0  while  the  server  idles;  thus  the  expected  change  in  v(B)  due  to  the  work-off 
is,  using  the  omni-convention  (i.e.  by  mentally  taking  the  expectation  of  each  side  of  the  equation): 

d^(B)  =  -dtP.ti'^u,)  (1.1) 

where  P.=  Pr(server  works)  =  fraction  of  time  as  seen  by  a  continuous  observer.  Of  course, 

P0  +  P,  =  l. 

(b)  New  arrivals  cause  upward  jumps  in  B.  A  pioneer  causes  a  jump  from  0  to  y,  and  a  follower 
causes  a  jump  from  w*  to  w,rx.  Thus,  Edi/>(B)  caused  by  new  arrivals  during  dt  is 

dt/>(B)  =  dt  AQo[^(y)-tl>(0)J  +  dt  AQ,[t/>(w,+x)-t/'(w»)]  (1.2) 

where  A  =  global  frequency  of  arrivals;  Q0=fraction  of  customers  who  find  the  server  idle;  and 
Q.=  fraction  of  customers  who  find  the  server  busy.  Note  that  the  Q- weights  refer  to  the  experience  of 
customers,  while  the  P-weights  refer  to  the  experience  of  an  observer.  In  a  model  with  a  steady 
Poissonian  source  of  customers  a  continuous  (or  random)  observer  sees,  statistically,  what  the  arrivals 
see  so  that  then  i>(u )  =  V'(w).  But  this  is  not  generally  the  case  for  a  non- Poissonian  source  (even  for 
a  piecewise  Poissonian  source  as  in  our  model)  where  customers  may  find  different  probabilities  than  do 
observers;  thus  tl’(u)  ^  0( w). 

From  (1-1)  and  (1.2)  we  get  the  balance  equation  for  the  backlog  in  G/G/l: 

*)  =  AQo[0(y)-il<(O)]  +  AQ*[i/’(w,+x.)-f’(w,)]  (1.3) 

Equation  (1.3)  can  be  integrated,  i.e.  brought  into  a  form  free  of  derivatives.  It  is  easy  to  show  (cf. 
Krakowski  September  1984  and  1985)  that 

U'(z)  -  tf(0)  =  z  tp'(z)  (1.4) 

for  any  positive  random  variable  z  which  is  interpreted  as  a  renewal  process;  z  is  the  residual  time  of  z, 
or  simply  its  residue.  Equation  (1.4)  plays  a  key  role  in  the  treatment  of  queues  when  integrating 
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omni-equations.  In  {act,  equation  (1.4)  helps  to  exploit  the  view  of  a  single- server  queue  as  a  renewal 
station  modified  by  gaps  in  service.  It  follows  from  (1.3)  and  (1.4)  that 


P.t/(u.)  =  AQ0y«/(y)  +  AQ.x«/(w.  +  x)  (1.5) 

Since  u>'(  )  is  also  a  general  function  of  its  argument  we  can  replace  in  (1.5)  each  ip1  by  v  —  we 

call  this  procedure  typographical  integration  —  thus  obtaining  for  G/G/l 

i/)(u.)=  AQ0yt(>(y)  4-  AQ.xi/>(w.+  x)  (1.6) 

From  (1.6)  we  get,  upon  setting  0(A)  =  1,  A(Q0y  +Q*x)  =  1  which  along  with  Q0+Q*  =  l  yields  Q0 
and  Q*  in  terms  of  A  and  y  and  x. 

Equation  (1.6)  can  be  also  read  in  the  following  manner:  The  random  variable  u*  is  a  mixture  of 
y  and  w*+x  with  weights  AQ0y  and  AQ*x.  Note  that  w*  and  x  are  independent,  as  follows 
from  the  derivation  of  (1.6);  unless  otherwise  indicated,  each  variable  in  each  argument  of  an  omni¬ 
equation  is  a  free  copy  of  its  generic  prototype:  it  is  distributed  like  this  prototype  and  is  independent 
of  any  other  variable  within  the  argument. 

The  unknown  processes  in  (1.6)  are  u.  and  w. .  The  variables  y  and  x,  and  hence  their  residues  y 
and  x,  are  known.  It  appears  that  the  theory  of  the  delay  in  G/G/l  and  in  some  of  its  variants 
|  reduces  to  relating  w*  to  u#  .  In  order  to  solve  (1.6)  for  both  w*  and  u*  we  need  another  equation  in 

these  variables,  wishfully  an  omni-equation.  However,  there  is  no  assurance  that  such  an  omni¬ 
equation  exists  for  all  or  most  models  with  general  input.  If  it  does  exist  its  form  is  likely  to  be 
complex.  (No  general  equation  relates  the  Laplace  Transforms  of  u.  and  w.,  which  equation  might  then 
be  transposed  into  an  equivalent  omni-equation.)  A  problem  of  applied  and  methodological  interest  is, 
for  what  models  can  w*  and  u*  be  related  by  means  of  omni-equations  alone? 

When  the  arrivals  are  Poissonian,  at  least  during  the  busy  period,  such  relations  between  u.  and 
w.  do  exist  and  can  be  simple.  From  these,  in  turn,  relations  between  u  and  w  can  be  derived. 
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Section  2.  Balancing  V'(B)  in  M/G/l  with  Exceptional  Service  and  Arrival  Rates 

Like  for  G/G/'l  we  find  it  useful  to  start  with  the  balance  of  v(B)  as  a  stepping  stone  towards  the 
analysis  of  the  delay  w.  Moreover,  the  backlog  is  of  economic  interest  when  some  costs,  e.g.  inventory 
costs,  depend  on  B. 

Let  f01  =frequency  of  transitions  of  system  size  from  N  =  0  to  N=l.  Clearly,  with  A=global  arrival 
rate,  and  with  a  and  A  being  the  Poissonian  intensities  during  the  server's  idle  and  busy  periods, 
we  have 


A  =  (tP0  +  AP.  =  AQ0  +  AQ.  ;  P0  +  P.  =  l 


From  (2.1a)  and  (2.1b)  we  get 


f0i  —  <rP0  —  AQ0 

AP.  =  AQ. 


and  Q0-fQ.  =  l 


(2.1a) 

(2.1b) 

(2-lc) 


Note  now  that  for  our  current  variant  of  M/G/l  we  have  the  key  property 

t(u.)  =  v)(w.)  (2.2) 

The  continuous  observer  who  sees  u«.  is  equivalent  to  a  random  observer  of  u,  who  takes  readings  with 
frequency  A  and  who  must  see  what  true  arrivals  see.  We  have  v( Poissonian  u.)  =  ^(continuous  u.) 
in  a  self-explaining  way. 


From  (1.6).  (2.1. b).  (2.1c)  and  (2.2)  we  get 

P.v(w.)  =  cryP0t(x0)  +  AxP.U'(w.  +  x)  (2.3) 

an  omni-equation  with  w.  as  its  only  unknkovn  random  variable.  Dividing  each  term  in  (2.3)  by 
P.  and  defining 

P  o  =  cry  and  p  —  Ax 


we  obtain 


t(w.)  -  (pa P0/P.)v(y)  +  ^(»'.+x 
With  ip(  )  =  1  in  (2.4)  we  find  that 

££-o=i-p 


X) 


(2.-1) 


P. 


and 


t’(w.)  =  (l  —  p)v’(y)-t-ptt'(w.  4-x) 


(2.6) 


Note  that  the  differentia]  equation  (1.3)  includes  in  its  arguments  y  and  x.  the  service  durations, 
whereas  the  integrated  equation  (2.6)  includes  in  their  stead  their  residues.  This  is  a  characteristic 
feature  of  integrating  a  differential  omni-equation  in  the  context  of  queueing  analysis:  Residual  service 
times  tend  to  go  with  integrated  equations.  This  explains  why  it  has  been  repeatedly  observed  that  the 
formal  structure  of  M/G/l  and  its  variants  is  simplified  when  expressed  in  terms  of  x  rather  than  x. 


From  (2.5)  and  Pn  +  P.  =  1  we  obtain 


Pn  = 


1  ~P 


[-P+  p  o 


and  P.  = 


and  from  (2.7)  and  (2.1)  we  find  A  and  Q0  and  Q.. 


1  —  P+  P  o 


We  find  the  successive  moments  of  w*  by  setting  ti-(A)  =  Ak  in  (2.6)  and  get 


(2.7) 


Ew.  =  Ey  + 


l-p 


Ex 


(2.8) 


With  c(A)  =  e  in  (2.6)  we  find  the  Laplace  Transform  of  w.: 


Ee 


-sw.  (1-P)E 


-sy 


1  —  pEe" 


(2.9) 


We  have  often  found  that  the  easiest  way  to  derive  the  Laplace  Transform  of  a  random  process  is  via 
an  omni-equation.  And  we  find  an  integral  equation  for  the  cumulative  distribution  function  of  w. 
by  setting  t  ( A  )  =  Pr{  A  <  t)  in  (2.6).  thus  obtaining  the  convolution  equation 
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Pr(w*<t)  =  (l-/>)Pr(y<t)  +  />Pr(w*+x<t)  (2.10) 

One  can  even  find  the  queue  size  by  specializing  the  function  or  functional  t>(A):  cf.  Section  5. 

The  reader  may  have  noticed  already  that  (1.6)  implies  (2.6)  when  V'(u,)  =  v(w*)  and 
\Q0  =  <tPq  and  AQ,=AP,  but  in  order  to  develop  the  M/G/l  model  independently  of  G/G/l  we 
started  again  with  the  basics. 
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Section  3.  A  Composition  Theorem  for  rp(w.) 

The  omni-equation  for  the  delay  w  in  a  regular  M/G/l  is.  in  integrated  form  (cf.  Krakowski  1986). 

ti’(w)  =  ( 1  —  p)V»(0)  +  pv>(w+x)  (3.1) 

R  R 

The  R  beneath  the  w  is  a  reminder  that  this  w  pertains  to  the  regular  M/G/l  . 

Theorem  1.  Suppose  that  certain  positive  random  variables  Z  and  g  satisfy  the  omni-equation 

V'(Z)  =  (1  -p)rp(g)  +  pip(Z+z)  (3.2) 

in  which  p  and  z  are  as  in  ( 3.1).  Then  the  following  composition  holds 

v(Z)  =  0(w  +  g)  (3.3) 

R 

which  says  that  Z  is  distributed  like  the  sum  of  the  generic  variables  w  and  g. 

R 

Proof  Shifting  (3.1)  by  g  we  obtain 


tf'(w+g)  =  (l-p)if(g)  +  pip(  w+x+g) 
R  R 


(3.4) 


From  (3.2)  and  (3.4)  we  get 


il>(Z)  —  pil'(Z  +  x)  =  t/>(w+g)  -  pV>(w -f x -f g) 
R  R 


(3-5) 


Defining 


<?(Z)  =  v(Z)  -  pi/^Z-fx) 


(3.6) 


we  g.  from  (3.5) 


Z)  =  d>(g+w) 

K 


(3.7) 


which,  <t>  being  a  general  function,  is  equivalent  to  (3.  3)  ;  this  completes  the  proof. 


Among  the  models  for  which  (3.3)  holds  is  the  multiple  vacation  model  in  whicvh  Z  =  w  and  g  = 

v.  a  vacation  period.  A  more  general  problem  than  the  simple  composition  in  terms  of  w  is  to  express 

R 

an  unknown  process,  e.g.  w  or  w*  or  B  or  B,,  in  terms  of  other  already  known  processes,  not 
perforce  the  delay  in  M/G/l  and  not  perforce  simple  convolutions. 
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From  (3.2),  (3.7)  and  (2.6)  we  have  the  composition 


V'(w.)  =  fc(w-t-y) 

_ _ K _ 


(3. 


Equation  (3.8)  solves  the  problem  of  deriving  the  delay  in  our  model  in  a  simple  and  versatile  wa\ 
We  are  now  ready  to  derive  the  omni-equations  for  the  delay  w  and  for  the  virtual  delay  u. 


1 


I 


i 


Section  4  The  Delay  in  M/G/l  with  Exceptional  Service  and  Frequency  of  Arrivals 


The  general  delay  w  is  a  mixture  of  two  delays,  pioneers'  and  followers',  with  weights  Q0  and  Q. 
as  perceived  by  the  source  of  customers.  The  delay  w  thus  clearly  satisfies  the  omni-identity 

i/'(w)  =  Q  0tl’(Q)  +  Q,ty(w.)  (4.1) 


Hence,  taking  account  of  (3.8),  i.e.  ip{ w.)  =  (w+y),  we  have 

R 


v(w)  =  Qo0(O)  +  Q.^w+y) 


(4.2 


which  is  a  valid  solution  for  w  since  all  coefficients  and  arguments  on  the  right-hand  side  of  the 
equation  are  known;  Q0  and  Q.  are  given  in  (2.1). 


The  steady  or  Poissonian  observer  will  write  an  omni-equation  for  the  virtual  delay  u  as  a 
mixture  of  pioneers’  and  followers'  components,  namely 

p(u)  =  Po^(0)  4-  P.^(u.)  (4.3) 

and  since  in  our  model  V>(w.)  =  t/>(u.)  we  can  also  write,  aided  by  (3.8), 


0(u)  =  P0^(0)  +  P.t/>(w  +y) 

_ R _ 


(4.4) 


Since  in  the  regular  M/G/l  ti’(u)  =  t/’(w)  we  can  also  write  (4.4)  as 

R  R 

tf’(u)  =  P0^(0)  +  P.t£>(u  +  y)  (4.4a) 

K 

It  is  seen  thus  that  the  omni-equations  for  w  and  for  u,  i.e.  equations  (4.2)  and  (4.4),  differ  only  in 
their  weights  for  finding  the  server  idle  or  busy. 


If  <7  =  A  then  P0=Q0  and  P.=Q.  and  we  have  ip(u)  =  t)’(w),  in  addition  to  ^(u.)  =  i/(w.), 
since  the  steady  (or  Poissonian)  observer  gets  the  same  statistical  picture  as  does  the  source  of 
customers,  both  pioneers  and  followers.  This  is  M/G/l  with  exceptional  service  but  uniform  arrival 
intensity,  dealt  with  by  Welch. 


In  some  other  variants  of  M/G/l  the  property  that  ^>(u.)  =  ti>(w.),  without  rp( u)  =  t/>(w),  has 


also  provided  the  key  to  the  treatment  of  the  queue  size  and  the  delay  under  steady-state  conditions. 
This  was  the  case,  for  example,  in  the  quorum  problem,  alias  Heyman's  A- Policy  (Cf.  Krakowski.  July 
1986  and  November  1986.)  Indeed,  when  analysing  a  variant  of  M/G/l  it  is  worth  considering  under 
what  conditions  observers  and  arrivals  see  the  same  statistical  picture. 


► 


Section  5  Poisson  Operator  and  Queue  Size 

Let  us  define  the  following  function  of  the  time  interval  z: 

#j(z)=  number  of  events  begot  by  a  Poisson  source  of  intensity  3  during  z. 

We  refer  to  as  the  Poisson  operator  of  intensity  3.  acting  on  z.  We  assume  that  the  Poisson  source 
is  independent  of  z  and  of  any  process  entering  the  model  unless  said  otherwise.  Where  no  confusion 
threatens  we  may  omit  the  subscript  3  and  the  qualifier  “of  intensity  3 ,”  thus  writing  #z  in  place  of 

#^(z)- 

The  distribution  of  #3(z)  is  well  known  (cf.  e.g.  Gross  and  Harris  1985): 


—  3z 

Pr(z=j)  =  E  - - jj - 

The  Poisson  operator  has  two  important  operational  properties. 

(1)  #(a  +  b)  =  #a  +  #b  if  the  intervals  a  and  b  are  independent,  in  particular  if  they 
overlap.  In  omni-notation 

v(#(a  +  b))  =  ^(#a+  #b) 


(5.1) 


do  not 


( 


(5.2) 


(2)  If  z  is  a  mixture  of  the  random  intervals  a  and  b  with  respective  weights  a  and  1— a  then  #z 
is  a  like  mixture  of  #a  and  #b.  In  omni-notation 


vi(z)  =  at/’(a)  +  (l  —  a)^(b)  =>  ^(#z)  =  a^(#a)  +  (l -»)^(#b)  (5.3) 

It  follows  from  the  above  two  properties  that  each  integral  omm- equation  with  constant  coeficients 
stays  valid  if  each  of  its  arguments  is  subjected  to  a  Poisson  operator  of  the  same  intensity.  We  can 
think  of  the  Poisson  operator  as  a  random  clock,  a  Poisson  clock,  which  assigns  a  random  discrete 
measure  to  any  time  interval  by  means  of  the  probabilities  (5.1).  < 

The  interpretation  of  the  resulting  omni-equation  has  to  proceed  most  carefully.  In  all  our 
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applications  the  Poisson  operator  will  have  the  intensity  of  the  customer  source  for  the  entire  process 
or  for  the  busy  period.  This  choice  often  leads  to  omni-equations  with  arguments  which  are  free  copies 
of  important  discrete  processes,  such  as  n  or  n,  . 

Of  course,  in  general  v(#(z, +z2))itty(#z1)-l- v(#z2),  even  if  z,  and  z2  do  not  overlap. 

Consider  now  the  omni-equation  for  the  regular  M/G/l,  i.e.  when  <7  =  A  and  t’(y)  =  and 

by  implication  Qo  =  Po< 

M/G/l  ti’(w)  =  (l-p)^(O)  +  p\p(w  +  k)  (5.4) 

Let  now  #=#A  and  let  ^(A)  =  o(#A)  for  any  argument  A  in  (5.4). 

Equation  (5.4)  becomes 

M/G/l  «H#w)  =  (l-p)<ji(#0)  +  d(#(w  +  x))  #  =  #A  (5.5) 

and  since  #(w  +  x)  =#w+#x  equation  (5.5)  becomes 

M/G/l  *(#w)  =  (l-p)*(#0)  +  <*(#w  +  #x)  #  =  #A  (5.6) 

It  is  clear  that  in  M/G/l  #w  has  the  same  distribution  as  qd ,  size  of  queue  left  behind  by  an 
entrant  into  the  service  station;  hence  ^(#w)  =  t£(qd) .  The  distribution  of  qd  is  like  that  of  n,  the 
number  of  customers  seen  by  a  random  observer  and  which  in  turn  is  distributed  like  qa ,  the  number 
of  customers  found  by  a  newcomer.  Thus  (cf.  Krakowski  1974) 

M/G/l  t£(#w)  =  t/j(qd)  =  v(qa)  =  tf’(n)  (5.7) 

Clearly,  #0  =  0  with  probability  1.  And  finally,  #x  =  number  of  Poisson  events  during  x  is 
assumed  known  along  with  x. 

We  therefore  find  that  (5.5)  becomes 
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M/G/l 


j  v’(n)  =  (l  — p)»'(0)  +  pv(n-t-h)  I  where  h  =  #x 


(5.8) 


The  power  of  the  Poisson  operator  to  recycle  omni-equations  for  the  delay  into  omni-equations  for 
the  queue  size  is  not  widely  known.  In  fact,  it  may  be  new. 


Note  It  can  He  shown  that  (cf.  Appendix)  with  #x  =  k  and  with  #x=h  we  have  the  omni¬ 
equation 


y.'(k)-y(0)  =  p[t/»(h  +  1 )  —  t/J(h )] 


(5.9) 


thus  expressing  the  functionals  of  h  in  term  of  the  functionals  of  k  which  have  been  more  often 
dealt  with  in  the  literature.  (Cf.  e.g.  the  k-matrix  in  Gross  and  Harris.) 


We  now  supplement  Theorem  1  of  Section  3  by  its  analogue  for  the  queue-size. 

Theorem  2  Suppose  that  the  non-negattve  and  integer- valued  random  variables  M 

and  ]  satisfy  the  omni-equation 

iHM)=(\-p)vU)+pri>(M+h)  h=#i  (5.10) 

m  which  p  and  h  are  as  tn  (5.8).  Then  the  following  composition  holds 


<l’(M)  =  4>(n+j)  (5.11) 

The  proof  is  analogous  to  the  proof  of  Theorem  1  and  shall  be  omitted. 

We  are  now  ready  to  recycle  (2.6),  the  omni-equation  for  the  positive  delay  w*  ,  into  an  omni- 
equation  for  n* ,  the  queue  size  while  the  server  works,  in  our  model  M/G/l  with  exceptional  service 
and  exceptional  arrival  rate. 

t’(w.)  =  (l—  p)v(y)  +  piJ}(Yi.+k)  (2.6)  =  (5.10) 

Applying  the  Poisson  operator  to  each  argument  in  (5.10)  we  obtain 

t(#w*)  =  (l-p)v(#y)  +  pV>(#w*  +  #x)  (5.11) 

Let  now  q„d=size  of  queue  left  behind  during  a  positive  delay;  q*a  =  £in  queue  found  by  a 
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follower;  n,=#in  queue  when  server  works. 

Clearly,  =  t/>(q,d);  =  ^(H*a)  follows,  as  can  be  shown,  from  the  fact  that 

transitions  “n=j  — *  n=j  +  l"  are  as  frequent  as  “n=j  +  l  —  n=j"  in  systems  where  all  events  are 
single  with  probability  one  (even  for  G/G/l  and  more  general  systems);  and  C(n,a)  =  t(n„)  because 
the  source  of  customers  and  a  continuous  observer  see  the  same  statistical  picture  in  a  queuing  system 
with  a  steady  Poisson  source,  and  in  our  model  the  source  is  steady  when  the  server  works.  (Cf. 
Krakowski  1973,  Theorem  B)  Therefore 


hum  (o.li)  and  (5.12)  we  obtain 


V>(n«)  =  (l  —  p)^(#v)  +  tfr(n«  +  h) 


From  Theorem  2  and  (5.13)  it  follows  that 


tl’(n*)  =  ti>(n+#y) 


We  can  also  derive  (5.14)  if  we  apply  the  Poisson  operator  to  each  argument  in  (3.8): 


^(w.)  =  tHg+y) 

R 


(5.12) 


(5.13) 


(5.14) 


(3.8)  =  (5.15) 


without  the  use  of  Theorem  2. 


We  have  clearly 


tf’fad)  =  Qo *H°)  +  Q*  v(q*d ) 

and  since  ti,(q,ct)  =  ti'(n*),  as  stated  in  (5.12),  it  follows  from  (5.14)  and  (5.16)  that 


<('(qd)=Qo^(0)+Q*^g+*y) 


(5.16) 


(5.17] 


Equation  (5.17)  applies  to  the  arrival  and  departure  streams  from  the  point  of  view  of  the  customer 
source.  Recall  that  t-(qd)  =  ti’(qa).  The  steady  observer  notes  that 
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v(n)  =  Pov-'(0)  +  P,v(n«) 


(•5.18) 


and  from  (5.14)  and  (5.18)  he  gets 


«'(n)=P0  ti'(0)  +  P,tf(n+  #y ) 

R _ 


(5.19) 


the  omni-equation  applicable  to  his  experience.  We  note  that  equations  (5.17)  and  (5.19)  differ  only  in 
their  weights  with  which  an  idling  or  a  working  server  is  found,  as  it  ought  to  be. 


page  19 


Section  6:  System  Size  and  System  Sojourn 

At  times  we  may  be  interested  in  the  system  size  and  the  system  sojourn  in  our  model  of  M/G/l 
modified  by  exceptional  service  and  arrival  rate  as  described  above.  We  can  easily  write  down  omni¬ 
equations  for  these  processes  from  the  results  for  the  queue  size  and  the  delay.  Thus  . 

^(Q)  =  Qo^(0)  +  Q.^(l+q*)=Qo«’(0)+Q*«-(l+n,)  (6.1) 

where  Q=system  size  found  by  a  newcomer.  Note  tnat  “system  size  tound  by  a  newcomer''  and 
“system  size  left  behind  by  a  served  customer”  have  the  same  distribution  (Cf.  Krakowski  1974). 

The  system  size  N  perceived  by  the  continuous  observer  is  described  by  the  omni-equation 

<*(N)  =  Potf(0)  +  P,iKl+n.)  (6.2) 

The  sojourn  time  W  in  the  system  is  described  by  the  omni-equation 

iA(W)=Q0tf'(y)  +  Q*r/>(w*+x)  (6.3) 

The  equations  (6.1)  and  (6.2)  and  (6.3)  are  self-explaining  as  mixtures  of  idle  and  busy  terms. 

Note  that  when  N>0  then  there  is  one  customer  in  the  service  station  and  n*  customers  in  the  queue: 
hence  the  argument  1  +  n*  in  (6.2).  A  similar  argument  holds  for  (6.1). 
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Appendix:  Relating  #x  to  #x 


Consider  the  omni-equation  of  the  renewal  process  generated  by  the  generic  service  time  x  (cf. 
Krakowski  1987): 

v(x)-r(0)=xv'(r)  r  =  x  (A.l) 

or 

w(x)  — ««(0)  =  lim  ^[v(r+dt)-v(r)]  (A. 2) 

Applying  the  Poisson  operator  to  each  argument  in  (A. 2)  we  have 

ti(#x)-V(#0)  =  lim  ^[v(#+  #t)-r/<(#)]  (A. 3) 

Denoting  #x  =  k  and  #r  =  #x  =h  we  write  (A.3)  as 


u(k)  — U-(0)  =  lim  ^[ti>(h  +  #dt)-tf’(h)j 


(A.4) 


It  is  clear  that 


9 

tp(h+ #dt)  =  Pr(one  Poisson  event)t/(h  +  l)  +  Pr(no  Poisson  event)\£>(h)  +  Odt 
=dtAr/)(h+ 1)  -f  ( 1  —  dtA)t/j(h)  +  Odt“ 


and  the  right-hand  side  of  (A.4)  becomes 

p(tt>(h-f  1)  —  t£(h)]  since  A x  =  p 


From  (A.4)  and  (A. 5)  we  have 

[~tr(k)-»’(0)  =  p[y.’(h+ 1)-  v(h)] 


(A.5) 


(5-9)  =  ( A.6) 
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